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ABSTRACT 

We use the recently released OGLE-III catalog of 17692 fundamental mode RR Lyr (RRab) stars 
in the Large Magellanic Cloud (LMC) to investigate the structure of its stellar halo. We apply 
conservative cuts in period, amplitude and magnitude to remove blends and other contamination. We 
use period-luminosity and period-color relations to determine distance and extinction of every star 
in our final sample of 9393 stars. In order to determine the scatter of our method, we compare the 
distributions of distances in two regions at the edges of the covered area with a central region. We 
determine the intrinsic line-of-sight dispersion in the center to be 0.135 mag or 3.21 kpc (FWHM of 
0.318 mag or 7.56 kpc), assuming zero depth in one of the edge regions. The conservative cuts we 
apply reduce the derived depth significantly. Furthermore, we find that the distribution of RRab stars 
is deformed in the sense that stars on the Eastern side are closer than on the Western side. We model 
the RRab distribution as a triaxial ellipsoid and determine its axes ratios to be 1:2.00:3.50 with the 
longest axis inclined by 6° from the line of sight. Another result of our analysis is an extinction map 
of the LMC and a map of internal reddening, which we make publicly available. 
Subject headings: galaxies: structure — Magellanic clouds — stars: variables: other 



1. INTRODUCTION 

The structure and evolution of the Large Mag- 
ellanic Cloud (LMC) stellar populations has been 
thorough ly investigated using many different rneth- 
ods fe.g. lvan der Marel et al.ll200l iNikolaev et al.ll200l 
lLah et al.ll2005t iCole et al.ll2005l ). One of the least stud- 
ied populations remains the stellar halo, for which the 
preferred tr acer objects are old m e tal-poor RR L y r vari - 
able stars (iKinman et alj 119911) . iMinniti et ail (|2003f ) 
and iBorissova et alj ( 2004 120061 ) measured radial ve- 
locity dispersions of a number of RR Lyr stars to be 
~ 50 km s" 1 , compared to e.g.. 20 km s _1 for carbon 
stars (|van der Marel et al.l 120021 ). This suggests that 
RR Lyr stars form a kinematically hot halo with a sig- 
nificant vertical spread. 

In this study, we use the OGLE-III ca talog of the fun- 
dame ntal mode RR Lyr (RRab) stars (fSoszvrisk i et al.l 
2009) to investigate the structure of the stellar halo of 
the LMC. The advantage of RRab stars is that distances 
to individual stars can be estimated using the period- 
luminosity relation and the extinction can be determined 
using the period-color relation. 

In Section [21 we describe the data set and how dis- 
tance and extinction are estimated for individual stars. 
In Section [3l we present the extinction map and the 
three-dimensional structure of the LMC. In Section [4] 
we discuss the possible systematics of our method, com- 
pare our findings to previous results on the structure of 
the LMC and discuss the implications of our results. 

2. DATA AND METHOD 

To study the structure of the LMC old stellar pop- 
ulation we use 9393 RRab stars. We start with 17692 
candidate RRab star s selected from the catalog of 
iSoszvhski et al.l (|2009f ) . who also provided the pulsation 
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mode classification based on Fourier coefficients of the 
light curves. We then exclude stars that were flagged as 
uncertain, foreground, or belonging to one of the LMC 
globular clusters. We applied further cuts to reduce 
the number of blended stars, likely misclassifications and 
stars with extreme properties. Specifically, we keep only 
stars with periods 0.45 days < P < 0.70 days, ampli- 
tudes 0.30 mag < AI < 0.85 mag and I > 18 mag. Fur- 
thermore, we remove the approximately 9% of stars that 
fall outside the main period-amplitud e sequence shown 
in Figure 3 of ISoszvhski et all (|2009f ). The equatorial 
positions (a, S) of the remaining 11171 stars are plotted 
in Figure Q] The edges of the OGLE-III coverage are 
rather irregular with respect to the RR Lyr distribution. 
To prevent appearance of artifacts in projected density 
plots we remove stars that fall outside a rectangle of size 
8.4° x 3.1° and oriented parallel to the long axis of RRab 
distribution on the sky. This leaves 9393 stars as our final 
sample. The rectangle is shown in black in Figure [Hand 
the method for determining the shape of the RRab dis- 
tribution is discussed in Section [3l The centroid position 
of our final sample is (ctQ,5o) = (80.35°, 69.68°). 

In order to estimate distances to individual RRab stars 
in the LMC we make use of the period-luminosity and 
period-color relations determined from RRab stars in the 
galactic globular cluster M3. We take advantage of the 
fact that the properties of RRab stars in this cluster are 
similar to the LMC, namely the mean period and the 
position of ridgelines in the period-amp litude diagram 
i|Alcock et aHl2000atlcTement et al.l l2008). Furthermore, 
the extinction toward M3 is very small, E( B — V) « 
0.01 mag (|Cacciari et al.ll2lM ISchlegel et al.lfl998h . For 
the slope of the perio d-luminosity relatio n 77 we use the 
value determined by IBenko et al.l ([2006), namely r\ = 
1.453. They also found that the scatter of their period- 
luminosity r e lation is 0.04 mag. From the dataset of 
IBenko et all (p006l ) we determined the relation between 
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Fig. 1. — Positions of the 11171 RRab stars satisfying cuts on period, amplitude and magnitude. Density contours marking levels of 250, 
500, 750 and 1000 stars per square degree are shown with grey lines. The rectangle containing the 9393 stars of our final sample is shown 
in black. Regions th at are investigated more closely are marked with blue, green and red parallel lines. The line of nodes measured by 
Nikolacv ct al. (2004) is shown with a black dashed line going through the centroid positi on of our sample (ctp, 8p). Black filled stars show 
positions of eclipsing b inaries with m easured distanc es; in order of incr easing a: HV 2274 (Ribas ct al. 2000; Grocncwcgcn & Salaris 2001), 
OGL E-0510 ifPiet rzvhski ct al. 2009), EROS 1044 URibas et al.ll2002T) , HV 982 HFitzpatrick et al J 12002) and HV 5936 ( Fitzpatric k et al.l 
120031) . 



period and intrinsic color (V — 7)o using robust linear 
least-squares fit. The period-color relation is 



[V-I) = 0.69 + 0.89 log P 



(1) 



with an rms scatter of 0.062 mag. 

The /-band extinction Aj to a particular star is es- 
timated from the difference between the observed color 
index (V — I) and the intrinsic color using 



Aj = R VI [(V - I) - (V - 1) ], 



(2) 



where Ryi = Ri/{R\r — Ri) is the ratio between total 
/-band extinction and reddening of the (V — /) color. 
If we get Ai < from Equation ([2|), we manually set 
Ai = 0. Throughout this paper, we set Ryi = 1.10, 
a slightly lower val ue than the usual Ryi — 1-55 (e.g. 
Schlcgel_el_aT|[l998), to avoid inconsistencies in extinc- 
tion estimates. We discuss implications of different val- 
ues of Ryi in Section [4TT1 

Using the information provided we are able to correct 
the observed intensity-mean / magnitude for extinction 
Ai and for period-luminosity relation according to 



I' Q =I-r,\ogP-Ai. 



(3) 



The line-of-sight distance z to a particular star is then 
obtained from /q assuming that the median distance of 
our final sa mple is 50 kpc . Our method resembles the 
approach of iNikolaev et al.l (|2004fl who studied the LMC 
Cepheids. 



3. THE THREE-DIMENSIONAL STRUCTURE OF THE LMC 

In this section, we first create an extinction map of 
the LMC and estimate the internal reddening. Then we 
focus on the three-dimensional structure of the stellar 
halo of the LMC as traced by RRab stars. 

3.1. Extinction map and internal reddening 

We construct the extinction map 1 of the LMC by di- 
viding the plane of the sky into pixels of 0.2° x 0.2° and 
calculating mean Ai for every pixel containing at least 
two stars. The extinction map is plotted in the upper 
p anel of Figure El Com parison with the extinction maps 
of Zaritsk v et alj (|2004f ) reveals the same basic structure. 
There is a high-extinction (Ai > 0.4 mag) region near 
30 Doradus extending toward north-west and a north- 
south filament of increased extinction (Ai ~ 0.25 mag) 
at a ~ 70°. Our map also shows the l ow ex tinction cen- 
tral region seen in the lZaritskv et al.l (2004) maps. 

Equipped with an estimated distance z to every star, 
we can estimate the line-of-sight structure of the extinc- 
tion. To this end, we select pixels with at least four stars 
and split the stars into two equal groups based on their 
line-of-sight distance. We calculate the mean Ai of each 
group and then subtract these two values to obtain the 
difference in extinction AAi between distant and close 

1 The extinction map can be obtained a t 
http : //www. astronomy . Ohio- st ate . edu/~ {}pej cha/lmc_extmap 
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Fig. 2. — Upper panel: Extinction map of the LMC. We calculated the mean value of Aj in pixels of size 0.2° X 0.2° containing at 
least two RRab stars. Otherwise, we assign Aj = 0. The blue, green and red regions from Figure [l] are marked here, too. Lower panel: 
Extinction variation along the line of sight AAj. We plot only pixels with A Aj > 0.05 mag. 
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Fig. 3. — Cumulative distributions of extinction-corrected RRab magnitudes I' for all 9393 stars (black solid line) and those in the regions 
marked in Figure [T] The Eastern region is shown with dashed blue line, central region with dotted green line and Western region with 
dash-dotted red line). The boxes near the center of the plot denote areas of 68% confidence areas around the mean. 



stars. The values of AAj are positive if the more dis- 
tant stars have higher extinction. The distribution of 
values of AAj is skewed to positive values with robust 
standard deviation of 0.063 mag and skewness of 1.48. 
We plot AAj as a function of position on the sky in the 
lower panel of Figure[2l As expected, the high-extinction 
regions seen in the upper panel of Figure [2] show up 
also in the internal reddening plot in the lower panel. 
A similar map of internal reddening was constructed by 
ISubramanian fc Subramaniaml (|2009f ) from the spread of 
colors of red clump stars measured by the Magellanic 
Cloud Photometric Survey. Their map also shows high 
internal reddening close to 30 Doradus, but the north- 
south high-extinction filament at a ~ 70° is not clearly 
visible. 

3.2. Structure of the Stellar Halo 

In order to determine scatter of our distance estima- 
tion method we first select two regions (the blue and red 
areas in Figure [1]) containing about 900 stars and lying 
near the edges of the covered area where we expect the 
intrinsic depth to be lowest. We can then attribute the 
measured width of the distance distribution to the inter- 



nal scatter of our distance estimation method. Cumula- 
tive distributions of I' for the selected regions are shown 
in Figure [3] with dashed blue and dash-dotted red lines, 
respectively. The widths of the distributions of I' mea- 
sured by the standard deviations are 0.108 mag (2.5 kpc 
at z — 50 kpc) and 0.123 mag (2.9 kpc) for the Eastern 
and Western regions, respectively, and are shown in Fig- 
ure [3] with blue dashed and red dash-dotted rectangles. 
We assume that the scatter of our distance estimation 
method is that of the Eastern region, i.e., 0.108 mag or 
2.5 kpc. We note that the scatter of period-luminosity 
and period-color relations gives 0.08 mag (1.9 kpc) as 
the scatter of our distance estimation method. Hence, 
our estimate of scatter in the distances is conservative. 

We expect the depth of the LMC to be highest close 
to the center. We plot in Figure [3] cumulative distri- 
bution (green dotted line) of I' of about 1100 stars in 
a 1° x 1° region centered on (ao,<$o)- The cumulative 
distribution is significantly wider than that of the side 
regions and also of the full sample that is shown with 
black line. The standard deviation of this subsample is 
0.173 mag (4.1 kpc) and is shown in Figure [3] with a 
green dotted rectangle. Using the scatter of our distance 
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Fig. 4. — Line-of-sight distance z of stars in our final sample as a function of rotated coordinates x' and y' . Blue lines mark contours of 
density 50, 100, 200, 300, 400 and 500 kpc" 2 for the (x',z) plane and 100, 200, 300, 400, 500, 600 and 700 kpc~ 2 for the (y',z) plane. Black 
stars mark the position s of the eclipsing binaries shown in Figure [T] note that we plot two distances for HV 2274 (Grocnewegcn fc Salarisl 
[20011: IRibas et aL 2000) and that HV 5936 falls outside the limits of the right panel. 



estimation method determined above, we can determine 
the intrinsic dispersion of depth in the central region to 
be 0.135 mag or 3.21 kpc (FWHM of 0.318 mag or 7.56 
kpc). As the intrinsic depth in the Eastern region most 
likely differs from zero, we can consider this value as a 
lower estimate. For comparison, the standard deviations 
of the distribution along axes of symmetry in the plane of 
the sky are 1.52 and 0.89 kpc, respectively. This means 
that the LMC RRab population depth is roughly twice 
the spatial extent in the plane of the sky. 

From Figure [3j we also see that the distributions for 
the Eastern and Western regions are shifted with respect 
to each other, in the sense that the Eastern region is on 
average about 0.05 mag (1.1 kpc) closer to us than the 
Western region. Indeed, stellar populations of the LMC 
are known to be a part of an disk with a line of nodes 
at th e position angle = 151.0° ± 2.4° (jNikolaev et all 
2004) and inclined so that the Eastern side is closer to the 
observer. From the RRab stars, we now see that the halo 
is also deformed so that the Eastern side is closer. En- 
couraged by these result, we now examine the full three- 
dimensional distribution of RRab stars. 

The population of RRab stars in the LMC can be ap- 



proximated by a triaxial ellipsoid. To determine the axis 
ratios and space orientation of the ellipsoid we first con- 
vert equatorial coordinates for every star to physical co- 
ordinates (a, S) — > (x, y) with the centroid of our sample 
corresponding to the coordinate origin, (ao, So) — > (0, 0), 
and keeping x to increase to the east. We then construct 
the covariance matrix of the coordinates (essentially the 
inertia tensor) and determine its eigenvalues and eigen- 
vectors. The ratios of the square roots of the eigenvalues 
are the ratios of the axes and the eigenvectors correspond 
to their spatial directions (for a similar analysis see e.g. 
iPaz et aLll2006|) . We first apply this method only in the 
plane of the sky described by the (x, y) coordinates and 
find that the distribution of RRab stars is elongated with 
an axis ratio of 1.71 and a major axis position angle of 
112.4°. We use this result to define new rotated coordi- 
nates, (x, y) — > (x',y'). Lines parallel to the longer and 
shorter sides of the rectangle in Figure[T]correspond to x' 
and y' axes, respectively. The results for the distribution 
in the plane of the sky are robust as they depend only 
on the precisely measured positions. 

Now we repeat the same procedure for the (x, y, z) co- 
ordinates, however, here the line-of-sight distance z is 
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Fig. 5. — Three-dimensional distribution of LMC RRab stars in 
(x',y',z) coordinates. The full distribution is displayed in semi- 
transparent gray and an isodensity surface of 56 kpc -3 is shown in 
dark yellow. The visualization was constructed from 0.5 X 0.5 X 0.5 
kpc~ 3 cubic voxels smoothed with a boxcar width of two elements. 

smeared by the internal scatter of our method. To correct 
for this, we subtract the value of internal scatter deter- 
mined above from the appropriate element of the covari- 
ance matrix. We obtain the axis ratios of 1 : 2.00 : 3.50, 
where the longest axis is inclined by 6° from the line 
of sight and the position angle of the projection of the 
ellipsoid on the plane of the sky is 113.4°. 

The finite sky coverage of our data can bias the deter- 
mination of axis ratios and directions. To see how seri- 
ous this problem is we remove stars lying outside the 250 
stars per square degree contour shown in Figure [1] and re- 
peat the procedure. We get the axis ratios of 1 : 1.36 : 3.53 
with the longest axis inclined by 3° from the line of sight. 
The decrease in the relative length of the second longest 
axis is caused by removing relatively more stars along x' 
than along y' . 

Figure [3] shows the spatial distribution of RR Lyr stars 
in the {x 1 , z) and (y', z) planes with overplotted density 
contours. Looking on the left panel we confirm our result 
that the Eastern part (with x' > 0) is on average closer 
to us than the Western part. In Figure [5] we give a more 
intuitive visualization of the spatial distribution of RRab 
stars in the LMC using a three-dimensional plot of an 
isodensity surface. 

4. DISCUSSION AND CONCLUSION 

In this section, we first discuss possible systematics 
that might be present in our analysis, then we compare 
our results to previous studies and conclude with impli- 
cations for the structure of the LMC and future work. 



4.1. Possible systematics 

There are a number of possible systematic errors that 
could lower the significance of our results. First of all, 
the assumption underlying all of our reasoning is that the 
population of RRab stars is homogeneous throughout the 
LMC. While we have no a priori reason to assume this, 
there are several checks for systematic errors. For ex- 
ample, comparing the period distribution of stars in our 
three regions, we find a mean period difference of about 
0.003 days, or about 0.01 mag in the period-luminosity 
relation. Another check is that the distributions of light 
curve Fourier coefficients show no differences between 
these regions above the level of 3% difference in the mean. 
In short, we do not see any evidence of different RR Lyr 
populations across the LMC. 

In our extinction correction of Equation l|2|). we set 
Rvi to a value lower than the fiducial Rvi = 1.55. 
We repeated our analysis with several values in the 
range 1.00 < Ryi < 1.55 and we found that the dis- 
tance distribution widths tend to increase with increas- 
ing \Rvi — 1 . 10 1 . For example, for Ryi — 1.55 we found 
standard deviations of 0.118, 0.118 and 0.191 mag for our 
Eastern, Western, and central regions, respectively. Fur- 
thermore, for values of Rvi > 1.2 we find inconsistencies, 
in the sense that stars closer to us have a higher aver- 
age extinction than the more distant stars. Or in other 
words, the mean AAi was found to be about —0.016 mag 
compared to 0.028 mag for Rvi — 1.10. Lower values of 
ratio of total to selective extinc tion to the LMC were 
found also bv lMisselt eFan (|1999f ). 

The cuts we applied on period, amplitude and mag- 
nitude tend to reduce the outliers and consequently the 
derived dispersion in depth. For example, applying no 
cuts yields standard deviation of 0.35 mag for the cen- 
tral region and 0.20 mag for the side regions; also the 
central region is on average significantly brighter than 
the side regions for all reasonable values of Rvi- Making 
our cuts even more stringent does not yield significantly 
smaller dispersions (or reduce the difference between the 
central and side regions). From this analysis we conclude 
that our cuts efficiently remove contaminated stars and 
that contamination of our sample is insignificant. 

We did not apply any cuts on color in order not to bias 
the derived shape in case there were large-scale extinc- 
tion variations. However, applying an additional color 
cut of 0.3 < {V — I) < 0.9 removes highly reddened 
stars, but does not change our results. Specifically, the 
intrinsic line-of-sight standard deviation in the central 
region would be 0.123 mag (2.91 kpc) and the axis ratio 
would be 1 : 1.99 : 3.14 with the longest axis inclined by 
9° from the line of sight. The application of the color cut 
changes our results on the intrinsic depth and axis ratios 
and directions only slightly. 

4.2. Comparison with Other Work 

iMinniti et aT] (|2003t) and iBorissova et all (|2004 120061) 
obtained radial velocities of 137 RR Lyr stars and found 
a very high velocity dispersion of about 50 km s _1 sug- 
gesting that the stars form a kinematically hot stellar 
halo. Our results on the depth dispersion are comparable 
with the results of lClementini et al. ( 2003) , who analyzed 
RR Lyr stars in two regions close to the bar and found a 
depth dispersion between 3.3 kpc and 3.8 kpc, especially 
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if we realize that the estimate of intrinsic scatter of our 
method is based on outer parts of the halo which prob - 
ably have substantial depth too. iSubramaniaml (2006) 
use d a smaller sample o f RR Lyr stars from OGLE- 
II (|Soszyhski et al.l I2003T ) to derive the scale-height of 
3.0 ± 0.9 kpc, which translates to dispersion of about 
6.0 ± 1.8 kpc. This is higher than our result, the differ- 
ence stems most likely from the fact that we estimated 
extinction correction on a star-by-star basis and we ap- 
p lied stringent cuts to r emove contaminations. 

iNikolaev et a l. (2004) found a number of Cepheids ly- 
ing significantly (> 7 kpc) both in front and behind the 
main disk of the LMC. Similarly, the distance determina- 
tions to eclipsing binaries i n the LMC show scatter larger 
than the uncertainties ( e .g.lGroenewegen fc Salarisll2001l : 
iFitzpatrick et al.l 12003 iPietrzvhski et al.l I2009T k Both 
types of objects are considerably younger and more metal 
rich than RR Lyr stars, but it is possible that the 
LMC halo contains a fraction of younger objects. A 
population of objects in a non-spherical halo elongated 
along the line of sight could help to explain the num- 
ber of observed mic r olensing events towards the LMC 
(lAlcock et al.l 120001* [Evans fc Kerinsl 120001: lAlvesI [20041 : 
iCalchi Novati et all 120061 : 1 Wvrzvkowski et al.ll2009h . As 
the microlensing cross-section depends basically linearly 
on the lens-source distance, doubling this distance with 
constant halo mass yields roughly twice cross-section of 
the self lensing due to halo of the LMC. 

4.3. Conclusion and future work 

We have investigated the structure of the LMC halo 
using fundamental mode RR Lyr stars. We found that 
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the halo line-of-sight depth is variable, in the sense that 
the depth is higher in the central areas than at the edges. 
Furthermore, the eastern side of the halo is in general 
closer to us than the western side, mimicking in part the 
LMC disk. We were also able to provide extinction maps 
of the LMC and resolve the internal extinction. 

While there is a significant increase of covered area 
between OGLE-II and OGLE-III catalogs, a survey of 
RR Lyr stars far from the axis of symmetry of their dis- 
tribution on the sky would help to better characterize 
their spatial distribution. It is possible that far from the 
center the distribution of RR Lyr stars might exhibit an 
overdensity similarly to Ceph eids, that are organi zed in 
a ring with an off-center b ar (Nikolacv et al. 2004|). 

Furthermore, iBorissova et al.l (|2006h found little vari- 
ation in the velocity dispersion as a function of distance 
from the center of rotation, but their sample has only 
two stars in t he farthest [2.0; 2.51 k pc bin. The OGLE- 
III catalogue ([Soszyhsk i et al.ll2009f ) gives stars as far as 
5 kpc from the center and it would be interesting to see 
whether the difference of the line-of-sight depth between 
the central and side regions appears also in radial veloc- 
ity measurements. 
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